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Abstract 

The thermal inflation is the most plausible mechanism that solves the cosmo- 
logical moduli problem naturally. We discuss relic abundance of superheavy 
particle X in the presence of the thermal inflation assuming that its lifetime 
is longer than the age of the universe, and show that the long-lived particle 
X of mass 10 12 -10 GeV may form a part of the dark matter in the present 
universe in a wide region of parameter space of the thermal inflation model. 
The superheavy dark matter of mass ~ 10 13 GeV may be interesting in par- 
ticular, since its decay may account for the observed ultra high-energy cosmic 
rays if the lifetime of the X particle is sufficiently long. 
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A large class of string theories predicts a number of flat directions, called moduli fields 
cj) |y]. They are expected to acquire their masses of order of the gravitino mass m 3 / 2 from 
some nonperturbative effects of supersymmetry (SUSY) breaking ||. The gravitino mass 
lies in a range of 10~ 2 keV-1 GeV for gauge-mediated SUSY breaking models || and in a 
range of 100 GeV-1 TeV for hidden-sector SUSY breaking models It is well known []5],[2| 
that such moduli fields are produced too much as coherent oscillations in the early universe 
and conflict with various cosmological observations. Therefore, we must invoke a mechanism 
such as late-time entropy production to dilute the moduli density substantially. 

The thermal inflation || is the most plausible mechanism to produce an enormous 
amount of entropy at the late time of the universe's evolution. In recent articles J!]] we 
have shown that the thermal inflation is very successful in solving the above cosmological 
moduli problem. Since it produces a tremendous amount of entropy to dilute the moduli 
density, abundances of any relic particles are substantially diluted simultaneously, which 
may provide a new possibility for a superheavy X particle to be a part of the dark matter in 
the universe .[] In this paper we show that it is indeed the case if the mass of X particle is of 
order of 10 12 -10 14 GeV and its lifetime is longer than the age of the universe. Such a long- 
lived X particle is particularly interesting, since its decay may naturally explain ||lT|-|T5H the 
observed ultra high-energy cosmic rays |]TB|-|T^] beyond the Greisen-Zatsepin-Kuzmin cutoff 



|19| when the lifetime is sufficiently long. 

In this paper we consider that the particle X was primordially in the thermal equilibrium] 2 ] 
and froze out at the cosmic temperature T = Tf when it was nonrelativistic Xf — mx/Tf > 1 
(i.e., the X is left as a cold relic. Then, the present relic abundance of X (the ratio of the 
present energy density of X to the critical density p cr ) is estimated by using the thermally- 



averaged annihilation cross section of the X, (ax \v\), as |22 



tt° h 2 = 0.76(72/ + l)x f 

X g,(T f ) 1 2M G (h^ Pcr /s )(a x \v\} 

where g*{Tf) ~ 200 counts the effective degrees of freedom at T = Tf, Mq ~ 2.4 x 10 18 
GeV denotes the reduced Planck scale, and so is the present entropy density (p cr /so) — 
3.6 x 10~ 9 h 2 GeV with the present Hubble parameter h in units of 100 km/sec-Mpc -1 . Here 
rif parametrizes the dependence on T of the annihilation cross section and rif — for s-wave 
annihilation, rif = 1 for p-wave annihilation, etc. Assuming the s-wave annihilation with 
(a x M) — m x 2 ( x f — 7 f° r m x ~ 10 13 GeV) we obtain 



1 Other candidates for the dark matter in the presence of the thermal inflation are the moduli 
themselves whose masses are less than about 100 keV f7|||], and the axion with a relatively high 
decay constant /„ ~ 10 15 -10 16 GeV gQ- 

2 For the superheavy particle X to be in the thermal equilibrium, the reheating temperature after 
the primordial inflation should be higher than mx ~ 10 13 GeV. Such a high reheating temperature 
is realized in e.g. hybrid inflation models [2C]. Although a large number of gravitinos are also 
produced in this case, they are sufficiently diluted by the thermal inflation and become harmless. 

3 Nonthermal productions of superheavy particles were discussed in Refs. |13,2l|. 



2 



tfxtf* . °^., =4.3xW-^) . (2) 



M G (ft- 2 ftr/.So) Vl0 13 GeV 

Therefore, the superheavy particle X of mx ^ 10 5 GeV leads to overclosure of the universe 
if its lifetime is longer than the age of the universe. In order to realize Vtxh 2 < 1 the 
dilution factor more than about 10 16 is required for mx ~ 10 13 GeV for example. However, 
such a huge dilution may be naturally provided by the thermal inflation. In this paper, we 
examine whether the thermal inflation could sufficiently reduce the energy density fix of the 
superheavy particle even if it was in the thermal equilibrium, as well as that of the string 
moduli. 

Let us start with reviewing briefly the thermal inflation model HPf. The potential of 
the inflaton field 5* is given by 

V = V -ml\S\ 2 + q^ r , (3) 

where -Wq denotes a soft SUSY breaking negative mass squared which is expected to be 
of order of the electroweak scale. M* denotes the cutoff scale and Vo is the vacuum energy. 
Then the vacuum expectation value of S is estimated as 

(S) =M= (n+ 2)"5Ph3J (m M™)^r, (4) 

and Vo is fixed as 

V = Ii±i„^, (5) 

so that the cosmological constant vanishes at the true vacuum. The inflaton a (^Re^ 1 ), 
which we call it a flaton, obtains a mass m\ = 2(n + l)rriQ. 

After the thermal inflation ends, the vacuum energy is transferred into the thermal bath 
through flaton decay and increases the entropy of the universe by a factor 



A ~ (6) 

Here Tr is the reheating temperature after the thermal inflation and determined by the 
flaton decay width r CT which can be written as 

r. = c.g, (7) 

where C a is a dimension- less parameter depending on the decay modes. We consider that 
the flaton o dominantly decays into two photons or two gluons, and C a is given byQ 



Here we have assumed that the flaton can not decay into two i?-axions which is an imaginary 
part of S to obtain a successful dilution. See Ref. @- 

5 We assume that the decay amplitudes are obtained from one-loop diagrams of some heavy 
particles. See also Ref. 0. 
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Then the reheating temperature is estimated as 

( 90 A 1 



7T 



r a M G ~ o.96 a 2 



1 mlM^ 



M 



(9) 



The entropy-production factor A, therefore, can be written as 

M 3 

A = s s r, 



130(n + 2)CfmlM| 



(10) 



and for n = 1 we obtain 

1.0 x 10 17 



A 
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Here the values m a = 100 GeV and M = 10 10 GeV correspond to M* = 3.5 x 10 18 GeV for 
n = 1. In the following analysis, we take n = 1 for simplicity. From Eq. (|TTD we see that 
the thermal inflation can dilute relic particle density extensively by producing an enormous 
entropy. Here it should be noted that the reheating temperature should be T R > 1 MeV to 
keep the big bang nucleosynthesis successful which leads to the upper bounds on M 

from Eq. (P) as 



M = 0.96 C3 



3 I 

1 ml M, 2 



G < 



T, 



R 



2.4 x 10 11 GeV 
3.9 x 10 12 GeV 



100 GeV 

m„ 



100 GeV 

These are translated into the upper bounds on as 

m„ 



for cr 



for cr 
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(12) 



AU = 3.2 < 



2.1 x 10 21 GeV 



5.4 x 10 23 GeV 



100 GeV 
100 GeV 



for a 
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(13) 



In the presence of the thermal inflation the relic abundance of the superheavy particle 
X [Eq. (H)] is reduced by the factor A in Eq. (|TT|) as 



6 In Ref. [E^] the lower bound on the reheating temperature is determined as about 0.5 MeV. 
Since our definition of the reheating temperature is different by a factor v3, it leads to Tr > 1 
MeV in our case. 
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5 

Q x h 2 = Q° x h 2 x I = 140 C| " >lx " > " 



A M 3 M£(h- 2 p cr /s ) 



°- 68 (ib^fev) 2 (lOb^ev) 2 (ib^kv) f ° r * ~> ^ 



(14) 



In Fig. |1] we show the contour of fix hi the parameter space of the thermal inflation model 
(in m a -M* plane). We find that the thermal inflation can naturally realize VLxh 2 < 1 keeping 
the constraint Tr ^ 1 MeV in a large region of the parameter space. 

As mentioned in the introduction, the thermal inflation was originally proposed as a 
solution to the cosmological moduli problem. Therefore, we turn to examine whether the 
thermal inflation could sufficiently dilute not only the density of the superheavy particle X 
but also that of the string moduli, simultaneously. 

When the Hubble parameter becomes comparable to the moduli masses, the moduli <p 
start coherent oscillations and the corresponding cosmic temperature is estimated as 

Ta ~ 7.2 x 10 6 GeV ( m ? T ) 3 . (15) 
VIOOkeV/ v ; 

Here notice that the moduli oscillations always begin after the X freezes out since Tf ~ 
m x 3> Tf/, even for the heavy moduli ~ m 3 / 2 — 100 GeV-1 TeV predicted in hidden 
sector SUSY breaking models. Because the initial amplitudes of the oscillations, 0o> are 
expected to be O ~ Mq, the present abundances^ of the moduli oscillations are given by 

Such a huge energy density of the moduli leads to cosmological difficulties for the typical 
moduli mass regions predicted in both gauge-mediated SUSY breaking and hidden-sector 
SUSY breaking scenarios. 

However, if the universe experienced the thermal inflation, the moduli abundances are 
reduced by the factor A [Eq. (0)] as @|] 

1 5 

2 \ mlrriaM, 




7 If the moduli masses are less than about 100 MeV, the moduli become stable until the 
present. On the other hand, for the heavier moduli mass region, $7^ is regarded as the ratio, 
i.P<t>l 's)d I '(Per I 'so)> where {p ( j ) /s)D denotes the ratio of the energy density of the moduli oscillations 
to the entropy density when the moduli decay. 
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We call these moduli produced at T = as "big-bang" moduli. In deriving Eqs. (|16D 
and ( |17|) we have assumed that the energy of the universe is radiation-dominated when 
the big-bang modulus starts to oscillate at H ~ m^. This assumption is justified for 
(m^/100 keV) 1//2 (mx/10 13 GeV) 2 ^ 2.9. On the other hand, when the energy at H ~ 
is dominated by the superheavy particle X, the present abundance of the big-bang modulus 



is related to the abundance of X (14) as 



m BB h 2 = flxh 2 [^-^ . (18) 

Therefore, both abundances are comparable for 0o ~ Mq- 

Furthermore, it should be noticed that the secondary oscillations of the moduli start 
just after the thermal inflation ends ||. We call the moduli produced by these secondary 
oscillations as "thermal inflation" moduli. The present abundances of the thermal inflation 
moduli are estimated as ||[7| 

mlM 



{%) TI h 2 = 6.0 x 1(T 2 C3- 



2.MZ(h-2n. \ M G 



mlM£(h-*p cr /s ) 
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L2 X ( lOOkevJ llOO^evJ llb^GevJ {W G ) ^ ° ^ 99 




• (19) 



We see from Eqs. (fP7|) and fll9 ) that the thermal inflation dilutes the moduli density 
substantially. In fact, it has been shown in Refs. PJ7|,f24]j that two moduli mass regions: (i) 

< 1 MeV and (ii) > 10 GeV survive various cosmological constraints in the presence 
of the thermal inflation. 

We are now at the point to examine whether the thermal inflation can solve the moduli 
problem and also dilute sufficiently the superheavy particle X at the same time. First, 
we consider the lighter allowed region of the moduli masses predicted in gauge-mediated 
SUSY breaking models. In Fig. ^| we show the contour plot of the abundance of X of mass 
rrix = 10 13 GeV as well as the various constraints in the m^-M* plane for = 100 keV. Such 
light moduli are constrained from the overclosure limit.0 Then the requirements (fl^BBh 2 ^ 
1 and (Q^Tih 2 < 1 put the lower and upper bounds on M* respectively. Furthermore, the 
condition Tr > 1 MeV leads to the upper bound on M* as represented in Eq. ([13|). We see 
from the figures that in the parameter space which survives above constraints one can obtain 
Qxh 2 ~ 10~ 6 -1 for mx = 10 13 GeV: i.e., the thermal inflation can dilute the abundance of 
the superheavy particle X sufficiently Similar results are obtained in the allowed moduli 
mass region ~ 10~ 2 keV-1 MeV. For example, when = 10 -2 keV, the thermal 



8 Notice that for the modulus of mass = 100 keV the constraint from the diffuse x(7)-ray 
backgrounds is not so severe as the overclosure limit. It gives a more stringent upper bound on 
the modulus abundance for 100 MeV ^ > 200 keV [25]. 
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inflation solves the moduli problem in the parameter regions m a ~ 10 2 GeV -1 GeV and 
M* ~ 10 13 GeV-10 17 GeV and in this region we obtain VL x h 2 ~ lO" 1 -! for m x = 10 13 GeV. 
We find that the thermal inflation is very successful, in a whole range of ~ 10 -2 keV-1 
MeV, to solve the moduli problem as well as to reduce the density of the superheavy particle 
whose mass is mx ~ 10 13 GeV even if it was in the thermal bath in the early universe. 

Next, we turn to discuss the case of the heavier moduli of masses m^ ~ 100 GeV-1 TeV 
predicted in hidden-sector SUSY breaking models. In Fig. |3| we show the contour plot of 
Vt x with various constraints for mx = 10 13 GeV and m^ = 100 GeV. Such heavy moduli 
are severely constrained not to destroy or overproduce light elements synthesized by the 
big bang nucleosynthesis. It has been found that (fl^BBh 2 and (fl^Tih 2 should be less 
than about 10 -5 ||26|| . Fig. [| shows that in order to solve the moduli problem, the thermal 
inflation requires M* higher than in the previous case, which results in more dilution of X, 
and hence we obtain VLxh 2 ~ 10 _9 -10 -7 for mx = 10 13 GeV. 

Although we have only considered the case n = 1, similar discussions also hold in the 
higher n case, except that the scale of M* becomes higher. 

In the present analysis we have taken the cutoff scale of the thermal inflation model 
M* as a free parameter. However, it is natural to choose it as the gravitational scale, i.e., 
M* ~ Mq. If it is the case, the thermal inflation dilutes the string moduli sufficiently only 
if their masses are m^ ~ 10 _1 keV -1 MeV. Moreover, as shown in Fig. f|, the dark matter 
density of the superheavy particle of mass mx — 10 12 -10 14 GeV becomes VLxh 2 ~ 10~ 4 -1, 
which is just the mass region required to explain the observed ultra high-energy cosmic rays. 
The required long lifetime of the X particle may be explained by discrete gauge symmetries 
or by compositeness of the X particle |28[ . 



In this paper we have shown that the thermal inflation provides indeed a new possibility 
that the superheavy X particle may form a part of the dark matter in the present universe. 
However, it gives rise to a new problem, since it also dilutes the primordial baryon asymmetry 
significantly. Because the reheating temperature of the thermal inflation should be quite 
low Tr ~ 1-10 MeV to dilute cf) and X sufficiently, the electroweak baryogenesis does not 
work. However, as shown in Ref. [29], the Affleck-Dine mechanism |3(J may produce enough 



baryon asymmetry even with tremendous entropy production due to the thermal inflation, 
if the moduli are light (m^ < 1 MeV) f25fl . Therefore, in the present scenario, the light 
moduli predicted in gauge-mediated SUSY breaking models are favored. 
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FIG. 1. The contour lines of the relic abundance (Qxh 2 ) of the superheavy particle X of mass 
mx = 10 13 GeV in the presence of the thermal inflation with n = 1 for the cases that a flaton a 
decays mainly into two gluons (A) and that a flaton a decays mainly into two photons (B). The 
contour lines are represented by the dotted lines. (The contour line of Qxh 2 =l is represented by 
the thick dotted line.) We show the upper bound on the cutoff scale M* obtained from Tr > 1 
MeV by the thick solid line. 
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FIG. 2. The contour lines of the relic abundance (Qxh 2 ) of the superheavy particle X of mass 
mx = 10 13 GeV in the presence of the thermal inflation with n = 1 for the cases that a flaton a 
decays mainly into two gluons (A) and that a flaton a decays mainly into two photons (B). The 
contour lines are represented by the dotted lines. (The contour line of Q,xh 2 =l is represented by 
the thick dotted line.) We show the upper bound on the cutoff scale obtained from Tr > 1 
MeV by the thick solid line. The lower bound on M* from (Q<h)BBh 2 < 1 and the upper bound on 
M* from (fl^Tih 2 < 1 for the case = 100 keV are also shown by the thick solid lines. 
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FIG. 3. The contour lines of the relic abundance (Qxh 2 ) of the superheavy particle X of mass 
nix = 10 13 GeV in the presence of the thermal inflation with n = 1 for the cases that a flaton a 
decays mainly into two gluons (A) and that a flaton a decays mainly into two photons (B). The 
contour lines are represented by the dotted lines. We show the upper bound on the cutoff scale M* 
obtained from Tr > 1 MeV by the thick solid line. The lower bound on M* from (f^ss/i 2 < 10 -5 
and the upper bound on M* from (Q^Tih 2 < 10~ 5 for the case = 100 GeV are also shown by 
the thick solid lines. 
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FIG. 4. The contour lines of the relic abundance (fiy/i 2 ) of the superheavy particle X in the 
presence of the thermal inflation with n = 1 for the cases that a flaton a decays mainly into two 
gluons (A) and that a flaton a decays mainly into two photons (B). We take the cutoff scale as 
M* = Mq ~ 2.4 x 10 18 GeV and = 100 keV. The contour lines are represented by the dotted 
lines. (The contour line of ^lxh 2 =l is represented by the thick dotted line.) We show the lower 
bound on the flaton mass m a > 1 GeV to open the gluon decay channel (A) and m a > 3.4 GeV 
from Tr > 1 MeV (B), and also the upper bound on m a from (Q^xih 2 < 1 by the thick solid 
lines. 
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